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We consider conformations and intra-molecular conformational transitions in amphiphilic star- 
like polymers formed by diblock copolymer arms with inner hydrophobic and outer polyelectrolyte 
blocks. A combination of an analytical mean-field theory with the assumption-free numerical self- 
consistent field (SCF) modeling approach is applied. It is demonstrated that unimolecular micelles 
with collapsed hydrophobic cores and swollen polyelectrolyte coronae are formed in dilute aqueous 
' solutions at high ionic strength or/and low degree of ionization of the outer hydrophilic block. An 

, intra-molecular conformational transition related to the unfolding of the hydrophobic core of the uni- 

molecular micelles can be triggered by a decrease in the ionic strength of the solution or/and increase 
in the degree of ionization of the coronal blocks. In the stars with large number of diblock copolymer 
arms the transition between conformations with collapsed or stretched core-forming blocks occurs 
continuously by progressive unfolding of the core domain. By contrast, in the stars with relatively 
small number of arms the continuous unfolding of the core is interrupted by an abrupt unravelling 
transition. A detailed SCF analysis indicates that under both unfolding scenario the arms of the 
star are extended fairly equally, i.e., no intra-molecular disproportionation occurs. 
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Q . I. INTRODUCTION 

c/3 ; 

Micelles of amphiphilic block copolymers attract a lot of attention in recent years because of their potential use in 

S applications including emulsion stabilization, viscosity regulation, biomedical applications, and nanotechnology 1 . One 
of the most interesting and promising fields of application of amphiphilic block copolymers is their use for controlled 
drug delivery. Many effective drug candidates have limited solubility, stability and can be even rather dangerous 
£Zj . because of their toxicity. To avoid these transport problems, drug delivery systems on the basis of amphiphilic block 
copolymer micelles with hydrophobic core and water-soluble corona have been proposed 2 ^. In this case the core 
encapsulates a drug whereas the hydrophilic corona makes this system soluble, protects from aggregation and even 
makes it selective with respect to the pore size thus allowing it to penetrate only into tumor tissues^. By attaching 
t— I \ a specific terminal group to the corona chain, biocompatibility and/or selectivity can be further enhanced. Upon 
^ ■ reaching the target compartment, a proper release of the drug molecules from the micelle should follow in response 
to the change in the surrounding conditions. The latter in most cases is the pH shift. For example, if the cellular 
uptake occurs via endocytosis 4 ^, the pH changes from its physiological value pH 7.4 to sufficiently lower value pH 
' 4-5 in lysosome. Drug release can occur via dissociation of the block copolymer micelle, cleavage of the corona chains 
or by conformational changes in the core itself. The use of starlike (AB) p block copolymers with hydrophobic inner 
B blocks and hydrophilic outer A blocks which form in solution the so called unimolecular micelles, Figure [1] a, has 
definite advantages over the self-assembled AB diblock copolymer micelles. Unimolecular micelles are more stable with 
I ■ respect to shear force, temperature, and pressure (often used for sterilization) than the block copolymers aggregate 
• • ' and, therefore, are more suitable for the noncovalent encapsulation of guest molecules. Moreover, formation of the 
. £h ' unimolecular micelles does not imply cooperative inter-molecular association, that is, there is no CMC. Therefore, 
vector systems based on the unimolecular micelles can be used at arbitrary low concentration. Another advantage 
is the prescribed "aggregation number", i.e. the number of star arms p and thus, a better control of the micelle 
size. This property is also very important for application of (AB) p unimolecular micelles as nanoreactors for synthesis 
of metal nanoparticles. At present, methods of synthesis of (AB) p stars are quite well developed^ and stars with 
different chemical nature of arm blocks solvable in water and organic solvents were obtained. They were successfully 
used for synthesis of goldi and platinum 8 nanoparticles. First steps in the creation of drug delivery systems on the 
basis of (AB)p unimolecular micelles were also made 4 . The final and, therefore, crucial stage in the successful work 
of a drug delivery system is the drug release when it reaches the target compartment. Drug release can occur, as it 
was mentioned above, via degradation of the core or of the micelle as a whole, but it may be useful to keep the (AB) P 
star intact (for example, if the core polymer itself cannot be safely removed from the body). This means that the 
drug release should be made through unfolding of the globular core, Figure [T]b, c. The driving force of this unfolding 
transition may be an increase of the repulsive interactions in the corona (or a decrease in attractive interactions in the 
core or both together) induced by variation of the environmental conditions. The corona of the (AB) p unimolecular 
micelle can be modeled as a convex spherical polymer brush grafted onto the B-core and immersed into a good solvent. 
This system was the subject of intensive theoretical studies in the last two decades beginning at the classical works 
of Daoud and Cotton^, Zhulin a 10 i 1:L , and Birshtein and Zhulina 12 . In more recent work of Zhulina, Birshtein, and 
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Borisov-i^ the Daoud-Cotton model have been revised, generalized, and extended on the convex polyelectrolyte brush. 
It should be noted that these approaches employed the so called "fixed ends" (or equal stretching) approximation, 
consisting in the assumption that the free ends of the grafted chains are localized at the edge of the brush, hence, all 
chains are equally, although not uniformly, stretched. For planar polymer brushes this approximation is often referred 
to as Alexander-de Gennes approximatio n 1 ^ 15 . 

The collapsed core of the unimolecular micelle can be also considered as a spherical polymer brush grafted to a small 
sphere with the free ends fixed at the outer boundary of the brush (hence, the core fits well to the Alexander-de Gennes 
approximation). Coil-globule transition in the star-shaped polymer in poor solvent was studied by Zhulina et al^ 
using scaling approach. It was demonstrated that under poor solvent conditions the intra-molecular concentration 
profile consists of two parts: a denser inner (central) zone is due to tethering of the star arms to the central point, 
here the polymer concentration decays inversely proportional to the distance from the star center. In the rest of the 
star the concentration is constant and is equal to the characteristic concentration of the collapsed polymer globule. 

The core-forming B-blocks of the amphiphilic star copolymer are subjected to the radial pulling force exerted 
by crowded coronal blocks. When this force exceeds a certain critical threshold, it may induce unraveling of the 
collapsed core and stretching of the B-blocks in the radial direction. In particular, the long-range Coulomb interactions 
between ionic (polyelectrolyte) A-blocks under low ionic strength conditions may be strong enough to provoke such 
kind of conformational transition. In 17 , conformational characteristics of amplhiphilic (AB)s star (with block lengths 
Na — 80, Nb = 20 and fraction of charged units in the A block Ub — 0.5) as functions of ionic strength were studied 
by means of numerical self-consistent field modeling. It was shown that at high salt concentrations all B-blocks are 
not extended by charged A-blocks and collapse in the core while at low salt concentrations the B-blocks are strongly 
stretched due to the corona pulling force. Monte Carlo simulation of (AB) p stars were performed by Nelson et. 
alr^ who considered not only the star with solvophobic core and solvophilic corona but the opposite interesting case 
(solvophilic inner + solvophobic outer blocks), paying more attention to the latter one. In these simulations, solvent 
was athermal for the solvophilic block, and the Flory parameter for the outer block was varied. It was found that 
the stars which have the innter solvophilic block can undergo a dramatic conformational change in which the outer 
solvophobic blocks aggregate into one or more compact globules. This aggregation transition is accompanied by a 
significant change in the size of the polymer as measured by the radius of gyration. In the case of stars in which the 
inner blocks are solvophobic and the outer blocks are solvophilic like in our study, no such transition was found. 

The aim of the present work, is to study equilibrium conformations of the amphiphilic (AB) p star block copoly- 
mers with hydrophobic B-blocks and strongly dissociating (quenched) polyelectrolyte A-blocks. We are particularly 
interested in the analysis of the stimuli-induced conformational (unfolding) transition in the (AB) P star triggered by 
the change in the ionic strength of the solution or/and degree of ionization of the A-blocks, which affect the strength 
of repulsive interactions in the corona. 

The rest of the paper is organized as follows: in Section "Model and general considerations" we specify our model of 
the (AB) p star block copolymer. In Section "Polyelectrolyte coronal domain" we calculate the pulling force exerted by 
the corona A-blocks onto the core blocks as a function of variable charge parameters and ionic strength of the solution. 
In Section "Core domain" we consider the B-core domain and analyze conformation of a collapsed in poor solvent star 
homopolymer subjected to radial extensional deformation. Our main results in terms of the diagrams of states and 
analysis of the character of the intra-molecular unfolding transition in the (AB) P star block copolymer are presented 
in Section "Conformational transition in (AB) p star block copolymer" . The analytical results are compared to the 
detailed self-consistent field (SCF) modeling in Section "SCF modeling" . In the end, we summarize our conclusions 
and present an outlook for the future work. 
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II. MODEL AND GENERAL CONSIDERATIONS 

We consider a dilute aqueous solution of (AB) P star copolymers. A star comprises p AB diblock copolymer arms, 
each consisting of a hydrophobic block B with the degree of polymerization Ng and a polyelectrolyte block A with 
the degree of polymerization Na- The arms are grafted by the ends of the B-blocks to a small sphere with the radius 
Tq. Both blocks are assumed to be intrinsically flexible. That is, the statistical segment length is of the order of a 
monomer unit length a. The block A comprises a fraction of permanently charged monomer units, i.e. is a quenched 
polyelectrolyte. Water is assumed to be a marginal good solvent for both charged and uncharged monomer units of 
the A-block. The short-range interactions between A monomers are modeled in terms of the virial expansion and 
only the pair monomer-monomer interactions with second virial coefficients va are considered to be relevant. On the 
contrary, for the hydrophobic block B water is a poor solvent, the short-range (van der Waals) interactions between 
B-monomers are modeled in terms of Flory-Huggins parameter xb, under poor solvent conditions \b > 0.5. Below we 
consider only very dilute solution of (AB) p star copolymer where the interactions between the stars are weak and do 
not affect their conformations. To develop an analytical theory we use a simple model assuming that all hydrophobic 
B-blocks are equally (but not uniformly) stretched so all the A-B junction points are located at the interface on the 
equal distance R from the center of the star. Then all hydrophobic B-blocks are within the sphere of the radius R 
while the coronal A-blocks form a swollen brush on the surface of this sphere. This approximation is verified and 
justified below by detailed assumption-free numerical SCF calculations. Equilibrium value of R corresponds to the 
minimum of the free energy of the star: 

F(R)=F A (R)+F B (R), (1) 

where Fx are contributions of hydrophilic (X = A) and hydrophobic (X — B) blocks. Minimization of the free 
energy, eq[l] is equivalent to the force balance condition 

f A (R)+j B (R)=0 (2) 

where fx(R) = dFx{R) / dR are the forces acting from the block X (X = A, B) on the phantom surface separating 
the blocks. In the absence of the corona hydrophobic B blocks collapse and form spherical core with the radius i?o- 
This core preserves its shape in the amphiphilic star if the corona pulling force fA(Ro) is not sufficiently strong. With 
increasing pulling force (for example, caused by a decrease in the ionic strength) it reaches a certain threshold value 
which is enough to surmount the core restoring force corresponding at least to partial unfolding of the hydrophobic 
core. Obviously, this threshold is higher for more hydrophobic core. As it will be demonstrated below, in a wide 
range of pulling forces, inner hydrophobic part of the (AB) p star has the shape of an "octopus" (compare with 19 ) 
with reduced dense spherical core and stretched B branches similar to extended hydrophilic corona branches. At very 
high pulling forces globular core disappears and the B-blocks acquire stretched in the radial direction conformation. 

In the following two sections we consider separately the swollen polyelectrolyte corona and hydrophobic core domain 
and derive the corresponding free energy as well as pulling and reaction force-radius dependences. The obtained 
dependences will be then used for analysis of equilibrium conformation of the (AB) p star as a whole. 

III. POLYELECTROLYTE CORONAL DOMAIN 

To calculate the corona free energy, we use the approach developed by Zhulina et alJ£. Consider a polymer brush 
composed by p chains uniformly grafted onto impenetrable sphere with the radius R. The free energy of the brush 
(per chain) can be expressed as 

F 3 f D fdr^ ' D 



k B T~2a 2 J R \dnJ dr + J R ^ mt ^ c P^ s ^ dr ' ( 3 ) 

where fi nt is the free energy density of interactions in k B T units. The latter equation implies that the end segments 
of all the A-blocks are located at distance D from the centre of the star, i.e. equal (but non-uniform !) stretching of 
the A-blocks is pre-assumed. The first term describes elastic stretching of the chain. The local chain extension dr/dn 
is related to the local polymer concentration c p (r) 

dn dn p 

p s(r) ■ dr dr 4nr 2 

s(r) is the area per chain at the distance r from the center of the grafting sphere: 

8 l r ) = — >r l . (5) 
p 
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The second term in eq[3]describes excluded volume and electrostatic interactions within the star and can be represented 
as 

fint {c p (r)} = fev {c p (r)} + f ion {c p (r)} (6) 

where 

fev {c p (r)} = v A c 2 p (r) (7) 

accounts for the excluded volume interactions in the corona. To specify the electrostatic contribution fi on {c p (r)} for 
ionic corona, the combination of the mean-field and the local electroneutrality approximation (LEA) is used. The LEA 
assumes that the local charge of polyelectrolyte coronal chains is compensated by the excess local concentration of 
counterions (including oppositely charged "own" polyelectrolyte counterions and salt ions). All ions in the system are 
monovalent. The LEA is applicable provided that the number of block copolymer arms in one star is sufficiently large 
so that the excess electrostatic potential is great enough to retain the mobile counterions inside the corona, even at low 
salt concentrations in the solution. In the LEA framework, the electrostatic interactions manifest themselves through 
the translational entropy of the ions, disproportionated between the interior of the corona and the bulk solution. 
The concentrations (the chemical potentials) of all mobile ions in the bulk of the solution are assumed to be fixed. 
The distribution of mobile co- and counterions between interior of the corona and the bulk solution is determined 
according to Donnan rule. The details of the approach can be found ir>i£. For strong (quenched) polyelectrolyte with 
constant degree of ionization a& it is given by-i£ 



fion {c p (r)} = - ^Jl + (a b c p (r)/§ ion ) 2 



a b c p (r) In ( \J 1 + (a b c p (r)/$ ion ) 2 - a & c p (r)/$ 



- 1 I $ion 

(8) 



where &i on = J2j c bj is the bulk concentration of all monovalent ions species. Therefore, the free energy can be 
expressed as 



F <- D 



R 



f{c p (r) ir } S (r)dr, (9) 



k B T 
where 

^• p > = ^-(£ a -^w + ^^> (10) 

is the free energy density. The outermost radius of the corona D is related to the polymer concentration profile c p (r) 
via the constraint of conservation of the total number of monomers 

D 4tt f D 

c p (r)s(r)dr = — / cJr) r 2 dr = N A . (11) 

R P Jr 

In the following we employ the so called local or quasi-planar (QP) approximation^ where instead of minimization 
of the functional, eq[9l with respect to c p (r) the local condition of vanishing the differential osmotic pressure in the 
corona: 

n (r)=c(r)——[ '-— =0 (12) 

5cp(r) \ c p (r) J 

is used. This leads to the following equation for c p (r): 



3 f P \ 2 1 , i<fi y/1 + (a b c p (r)/<S> ton )2 - 1 



from what the polymer density profile follows. Explicitly it can be found only in the "inverse form" r = r(c p ) 

1 /*<% A [V A <% + $ionWl + (a bCp /$ ion ) 2 - 1)]V 



r 



(JL) 1/2 31/4 (14) 
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FIG. 2: Free energy of polyelectrolyte corona (per chain) (a) and corresponding pulling force (b) as functions of the core radius 
for different salt concentrations, cj>i 0n = <I>; on a 3 . The values of parameters are Na = 1000, p = 20, at, = 0.5, va = a 3 /6 



Then, the normalization condition, eg 1111 can be re- written as an integral over c p : 



47T 



D 



47T 



— / c p (r) r 2 dr = — I c p r 2 (c p ) dc p = N A 

P Jr P Jc p (R) dc p 



The free energy functional can be written as 



F A 



-in 



f{ c p( r )> r } r2 dr 



4tt 
P 



k B T P Jr 

where f(c p ) is obtained from / {c p (r),r}, eqUHl by using eg [HI 



Cp{D) 



f(c P y( Cp )^ldc p 



(15) 



(16) 



f(Cp) 



1 



-VACp 



2 2 
— ctbC p In 



yjl + (a b c p /<$> ion ) 2 - ] 

ton -I o n J 



(17) 



We can obtain the dependence of the free energy Fa as a function of the core radius R as follows. Polymer concentration 
at the core surface, c p (R), is related to the core radius R via eqQ2]or cqQJJ Once parametrically c p (R) is known, 
the polymer concentration at the periphery of the brush, c p (D), is found using normalization condition, eq 1151 then, 
finally we use egs \W\ and [171 to find Fa- The stretching (pulling) force from the corona side applied to the B-blocks 
at the A-B junction points is expressed as 



u 



dF A 
dR 



kuT- 



4-7T 
P 



f{ Cp (D)}D 2 — -f{c p (R)}R 2 



A<jr 

-k B T-—R 2 f{c p (R)} 
P 



f{c P (D)} c p {R) 
f{c p (R)} c p (D) 



(18) 



dD/dR is found by differentiating the normalization condition eq [TT] with respect to R. Figure [2] shows the free 
energy and pulling force dependences of the radius R calculated for Na — 1000, p — 20, a b — 0.5, va = a 3 /6 and 
various volume fractions of salt, 4>i on = $i on a 3 . Both Fa and are decreasing functions of R, a decrease in the ionic 
strength leads to an increase of the pulling force jU. 



IV. CORE DOMAIN 



The core domain of the (AB) p star can be considered as a polymer star under poor solvent conditions with each arm 
subjected to radial tension exerted by the swollen coronal blocks. The analysis of the conformations of the core blocks 
necessitates taking into account the specific effects appearing upon extensional deformation of a polymer globule. The 
latter reveals, in its turn, an non-trivial physical picture. 
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A. Extensional deformation of a single chain in poor solvent 

As it was first shown by Halperin and Zhulin a 20 ' 21 , when a constant force is applied to the ends of a polymer chain 
collapsed in poor solvent, there exists a critical force, below which the globule remains compact but changes its shape 
from spherical to prolate ellipsoidal, whereas above the critical force it unwinds into a stretched chain. When the 
globule is deformed in the extension ensemble, i.e. the distance between chain ends is imposed, the unfolding scenario 
is more complex: it occurs via formation of the "tadpole" conformation in which the spherical globular "head" coexists 
with the stretched string of thermal blobs ( "leg" ). In our recent work o 22 ' 23 we have performed a detailed self-consistent 
field (SCF) calculations and developed a quantitative analytical theory of the equilibrium unfolding of a globule formed 
by a flexible homopolymer chain collapsed in a poor solvent and subjected to an extensional deformation. Both the 
SCF calculations and the analytical theory prove that upon an increase in the end-to-end distance three regimes of 
deformation successively occur: at small deformation the globule acquires a prolate shape, the reaction force grows 
linearly with the deformation; at moderate deformations the globule is in the tadpole conformation, the tail of the 
tadpole is uniformly stretched, its extension as well as the reaction force is weakly decreasing with the deformation; 
then at certain extension the globular head unravels and the globule completely unfolds, the reaction force drops down 
and then grows again upon further extension (this jumpwise "unraveling transition" was first predicted by Cooke and 
WilliamsSi). In the framework of the developed theory it was possible to calculate prolate globule-tadpole and 
tadpole-open chain transition points and to find corresponding conformational characteristics, such as the asymmetry 
and the number of monomers in the globular head, reaction forces and force jump at transition, in a wide range of 
polymerization degree N and solvent quality (%). Our analysis has shown that the system exhibits a critical point; i.e., 
there exists a minimal chain length N cr (x) below which, at N < N cr (x)i the intramolecular microphase segregation 
in the extended globule does not occur. The globule is deformed as a whole, without intramolecular segregation but 
by progressive depletion of its core. 



B. Homopolymer star in poor solvent 

Consider a free polymer star comprising p arms each of Ng monomer units immersed into a poor solvent, \ > 0.5. 
The star can be represented as a spherical polymer brush grafted onto a small sphere with the radius ro- The mean 
grafting area per chain is then sq = Air^/p. Since p arms cover the surface of the sphere and the arm cross-section 
area equals to a 2 , r can be estimated as r = ayp] (47r). It is assumed that all the chains are extended equally but 
non-uniformly in the radial direction so that their free ends are localized at the edge of the star at the distance Rq from 
its center. It can be showni^ that the concentration profile of the collapsed star consists of three concentric regions: 
Close to the center, there is a narrow zone where polymer volume fraction is approximately equal to unity. Then it 
decays as 1/r down to the concentration tp which is the characteristic globular concentration. In the Flory-Huggins 
model, the polymer concentration tp or, equivalently, the volume fraction in the globule (f> = tpa 3 can be expressed 
via Flory parameter \ as ^ 3 " X — ~ ^^T^ — ^ ■ Since the radius of the central zone of decaying concentration is small 
compared to the radius of the globule Rq, in the following calculation we assume that the polymer concentration in 
the collapsed star is uniform and equal to tp, i.e. 

This simplification should not affect, neither qualitatively nor quantitatively, final result. Consider now the situation 
where the ends of star arms are "pinned" at the distance R > R , i.e. each arm of the collapsed star is radially 
stretched. In accordance with the general picture of polymer globule deformatio n 20 i 22 we propose the following model: 
we assume that under such extensional deformation, the star segregates into corona formed by equally stretched " tails" 
and spherical globular core, Figure [3] b. The polymer concentration in the core is still equal to tp whereas the tails are 
the uniformly stretched parts of the polymer chains exposed to the solvent. If the core contains N g i monomer units 
per arm and has, therefore, the radius 



3N glP \ 



1/3 



then the free energy of the deformed star per star arm can be represented as follows 

ft. _ A r (<t\ ir+llNil+1 .^ + w - Ml (21) 



k B T 2a 2 J ro \dnj ^ 91 ' p 2(N B ~ N gl )a 2 



(a) (b) (c) 

FIG. 3: Polymer star collapsed in poor solvent: (a) free (unperturbed) star, (b) radially stretched star at moderate deformation, 
(c) strongly radially stretched star. 



where the first term is the elastic free energy of the globular core (r is the radial coordinate and n is the monomer 
ranking number counted within each arm from the center of the star), k B T/i is the monomer chemical potential 
(free energy) in the infinite globule (fi is negative under poor solvent conditions), the third term accounts for the 
intcrfacial contribution from the surface of the globule, k B T^ is the interfacial tension coefficient (7 is positive due 
to the energetic and conformational penalty of redundant monomer-solvent contacts at the interface), and the last 
term is the elastic free energy of the tail. In the volume approximation, which is applicable for sufficiently large core, 
it is assumed that both /1 and 7 are functions of the solvent quality only (i.e. independent of the core size), their 
dependences on the Flory parameter x are expressed as follows: the monomer chemical potential \x can be found, 
once the polymer volume fraction within the globule, 4>, is known^ 3 -: 

/i = 2 + ^log(l-£). (22) 

For the interfacial tension coefficient k B T-f an explicit expression can be obtained for moderately poor solvent onl y 23 ' 25 : 

7« 2 -^(l-2x) 2 (23) 

Comparison with the values of 7 obtained numerically using Scheutjens-Freer self-consistent field (SF-SCF) approach 
shown^ 3 - that this approximation works good for x < 1.2 which is the range of \ we are interesting in. Since the 
polymer concentration in the globule is constant and equal to ip, local chain tension is dr/dn = p/(4irr 2 ip), and the 
integral in eq[5T]can be easily calculated 

FB _ 3P fin +T ,*a + »<*-y' ,24, 



k B T 8na 2 p\r Rgi J 4 ' V 2(N B - N gl )a 

To find the equilibrium value of N g i at given R, the free energy, eg 1241 should be minimized with respect to N g f. 

dF B _ dF B 8F B dR gl _ 

dN gl dN gl dRgi ' dN gl ' 1 ' 

which leads to the following equation 

3P 2 2 MR-R g i) 3(i? - R g if 

32^Vi?4 ;fl 2 +7 vRgl A^R 2 gl (N B -N gl )a 2 + ^ ^ 2{N B - N gl fa 2 ' 1 ' 

Solution of eq [55] should give the equilibrium number of monomer units (in a star arm), N^, that remain in the 
globular core as function of the imposed star deformation R. Eq [26] can be resolved analytically with respect to 
R = R(Ngi). The solution is as follows: 



R = R g i + (N B - N g i)a 



p 2 47 




47T(y5_R~ i a y 3 3ipR g i 



(27) 



where R g i = R g i (N g i ) is given by eq [501 Let us consider eqs [27] and [5D] as an implicit dependence of the number of 
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FIG. 4: Number of monomers in the star core N g i as function of the star deformation R calculated using eq[27jfor \ — 1, 
Nb = 200 and various values of p — 10, 20, 30, . . . , 80 in the main figure and p = 10, 12, 14, . . . , 20 in the inset. Dashed parts 
of the curves correspond to unphysical solution. 



monomer units in the core, N g i, on the given value of the star extension R. Figure [4] shows a family of N g i vs R 
curves for Ng = 200, x = lj an< i various number of arms p. Solid lines show solutions corresponding to N g i that 
decreases with increasing extension R whereas dashed lines show unphysical situation (simultaneous increase of the 
deformation R and of the amount of the polymer condensed in the globular core) and should be therefore excluded 
from our consideration. The point where dR/dN g i = represents a spinodal point where the two-phase "octopus" 
state (Figure [3] b) looses its stability. It corresponds to the (i) largest possible extension and (ii) smallest possible 
core size for the locally stable microsegregated deformed star. It can bee seen that for moderate values of p the 
smallest number of monomer units in the core N g i m in > 1, for instance at p = 20, N g i t mm ~ 2.7. As p grows N g i t m in 
decreases and assumes values less than unity. 

By substituting eq [37] into eq one obtains the equilibrium free energy of the segregated deformed star (Figure 
[3] b) as function of N g i ; the latter expression together with eq [37] provides a parametric dependence of Fb (R) ■ This 
free energy should be compared with that of completely unfolded star which has no globular core (Figure [3] c), 



F, 



[unfolded] 



B 



HR - r f 
2N B a 2 



(28) 



to decide which state of the star (partly or completely unfolded) is thermodynamically stable. The comparison of the 
free energies of the star with partially and completely unfolded arms enables us to localize the first-order transition 
point between these two conformations and to calculate the magnitude of the jump-down in the reaction force. As 
follows from results presented in Figure an increase in the number of arms leads to the shift of the transition point 
towards larger values of deformation R and to a decrease in the magnitude of the force jump. The reaction force fs 
per arm of the partially unfolded core is: 



f B (R) 



dF F , 



dF B dF B ON, 



gl 



dR dR dN g i dR 
Using eq[27]we can express fs as a function of R g f. 



3(R - R g i) 
{N B ~N gl )a2 



fe(Rgl) 



4mpR? gl a 



47 

3(pR gl 



The reaction force per arm of completely unfolded/unraveled star is 



Aun folded) 
IB 



(R) 



dF 
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B 

dR 



HR-rp) 
N B a 2 



(29) 



(30) 



(31) 
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FIG. 5: Equilibrium force-extension curves for \ = 1> Nb = 200 and various values of p — 10, 20, 30, . . . , 80 in the main figure 
and p = 10, 14, 18, . . . , 30 in the inset. 



Figure [5] shows equilibrium force-extension curves calculated for the set of parameters (N — 200, x = 1) similar to 
that in Figure 2) It can be seen that if the number of star arms is small (p = 10 or 14 in Figure [5]), the force-extension 
curve looks very similar to that for a single chain 22 ' 23 . There is, however, a tiny difference in the initial part of the 
force-deformation curves: For the star the curve starts at R = Ro with a non-zero force. This is because in the present 
work we do not consider an initial deformation regime that could lead to a distortion of the spherical symmetry of 
the collapsed star or to a decrease in the mean polymer concentration (p within the star. Hence, a pulling force which 
is less than the threshold force that gives rise to the microphase segregation does not affect the structure of collapsed 
star. The value of this threshold force can be easily obtained from eqEOlby setting N g i = Ng and R g i = Rq : 



P 



A-KLpR^a 



47 



(32) 



One can see that the threshold force grows with increasing number of arms in the star p. Similar effect is observed 
for a globule formed by a linear macromolecule: the value of the force corresponding the onset of the microphase 
segregation in the deformed globule grows as the degree of polymerization increases. In the star case the overall 
molecular weight is Np and this causes the effect. In the phase segregation regime, the force for the star with p = 10 
weakly decreases as a function of the arm extension, in the "unraveling" transition point it drops down and then 
grows as 3(i? — r^)/ (Nbcl 2 ) ~ 3R/(Nbo. 2 ) for the unfolded star. However, with an increase in p the character of the 
force-extension curves changes: the slope of the quasi-plateau on the fs(-R) curve becomes less negative and can even 
be positive. To estimate the slope of the force-extension curve, the derivative dfs/dR should be calculated. Since the 
reaction force is given in eq [301 as a function of R„i , it is expressed as follows 



dfs = 9fs_ dR^ dN^ 
dR dR g t ' dN g i ' dR 

Since dR g i/dN g i > and dN g i/dR < 0, the sign of dfs/dR is opposite to that of dfs/dRg 
can be easily found: 



(33) 



Expression for dfs/dRgi 



d 



dR g i \k B T ) 



2-KipR 



3Ki 



20 



47 



(34) 



We see that at small p the second positive term in square brackets in eq [M] related to the globule's surface tension 
dominates, then dfs/dR < 0, and the quasi-plateau has a negative slope, whereas for large p the first negative 
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FIG. 6: Free energy of (AB) P core-shell star with p = 20, Na = 1000, iV s = 200, v A = a 3 /6, \b = 1, a b = 0.5 and different 
values of <j>ion- Arrows show the position of local or boundary minima, dashed line separates regions of partially and completely 
unfolded core. 



term dominates and this corresponds to the increasing quasi-plateau. The threshold value of p = p* at which the 
negative-to-positive change in the quasi-plateau slope occurs can be easily found by considering the initial part of the 
force-extension curve where N g i — Nb, R g i = Ro- H P = P* corresponds to dfs/dRgi = 0, this leads to the following 
equation 

P , 2 7 



(35) 



We can find an explicit solution of the latter equation if we take Rq ~ (^Nsp/^ip) 1 ^ 3 and |/x — 3^r — \J ~ f M- 
With these simplifications we obtain 

p=47r U^J F^F' (36) 

Hence, p* increases upon an increase in the solvent strength for the B-block. For the case Nb — 200, \ = 1 presented 
in Figure U] and Figure [SJ eq [321 gives p* ~ 25 which is in accordance with the results presented in the inset in 
Figure [5] Note also that according to eq[36]p* scales as a square root of arm length Nb- At the same time, as p 
increases, the drop in the reaction force that accompanies the unraveling transition decreases. This is in accordance 
with the conclusion made for a deformed single globule where the force jump at the tadpole - uniformly stretched 
chain transition decreases as the overall chain length N increases. 



V. CONFORMATIONAL TRANSITION IN (AB) P STAR BLOCK COPOLYMER 

In the preceding sections the expression for the core and the corona free energies and corresponding forces were 
obtained. Consider now the (AB) P core-shell star as a whole. Figure [5] presents as an example the free energy 
of the star F(R) = Fa(R) + Fb(R) as a function of the radial position R of the A-B junction points at different 
salt concentrations (Fa(R) and Fb(R) were calculated using expressions obtained in the previous sections). An 
equilibrium conformation of the star corresponds to a (global) free energy minimum. The latter is equivalent to the 
core and corona force balance condition, eq[21 i.e. alternatively we can consider force f vs radius R dependences. Such 
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R/a 

FIG. 7: Core reaction force (gray curve) and corona pulling force (black curves) in (AB) P core-shell star withp = 20, Na = 1000, 
Nb = 200, va = a 3 /6, \b = 1, <*b = 0.5 and different values of 4>ion- Circles show the position of equilibrium between the core 
and the corona. 

dependences corresponding to the case presented in Figure |H] are shown in Figure [JJ for the core, there is a single 
curve starting at r = Rq while for the corona, there is a family of curves corresponds to different salt concentrations. 
Figure [JJ can also be considered as a (partial) superposition of Figure [2] b and Figure [5l 

At relatively high salt concentration, <pi on = 10 _1 the electrostatic repulsions are weak, and, as a result, the pulling 
force from the corona with the inner radius i?o, eq[TJJl is smaller than the threshold force, eqEU necessary for unfolding 
of the core, We see that in this regime the free energy has one boundary minimum corresponding to a unimolecular 
micelle with collapsed core (R = Ro = (ZNp/l-Ky) 1 /*). Correspondingly, the \f A (R)\ curve does not cross fs(-R) 
dependence, as it can be seen in Figure [JJ 

At lower salt concentration 4>%on = 10 -2 , which is below certain salt concentration threshold, |fyi(-R)| curve crosses 
fs(-R) curve in the microphase segregation ("octopus") regime at R > Rq. This point corresponds to the local (non- 
boundary) minimum of the free energy F(R). A similar picture is observed at <fii on = 10~ 3 , the equilibrium radius 
of the B-domain, R, is shifted towards larger values. At 4>i on — 10 the free energy minimum and the point of 
intersection of and fs(-R) curves correspond to the conformation with completely unfolded core. Remarkably, 

close to the transition point between the conformations with partially and fully unfolded core, the free energy has 
two minima (this can be seen on curves for 4>i on = 10 -3 in Figure [6j in Figure [7J two points of intersection are clearly 
seen). Besides the main minimum corresponding to the octopus- like conformation, the second metastable minimum 
related to completely unfolded core appears. The appearance of the two minima in the free energy implies co-existence 
of the stars with the partially and completely unfolded core domain and signifies quasi-first order character of the 
unfolding transition: Upon a decrease in 4>i on the system "jumps" from the "left-hand-side" minimum (corresponding 
to the microphase-segregated partially unfolded core) to the "right-hand side" minimum (corresponding to the star 
with completely unfolded core domain). For the set of parameters used here, the co-existence between microphasely- 
segregated and fully unfolded core occurs at R s=s 61, as it is shown in by a vertical dashed line in Figure [6] Figure 
[5] shows equilibrium position R of the A-B junction points, the overall star radius D (which is, equivalently, the 
outermost corona radius), and the corona thickness H = D — R as functions of salt concentration <f>i on - As we have 
shown, with a decrease in <pi on the star undergoes a sequence of intramolecular conformational transitions: when 
salt concentration is large, repulsive interactions in the corona are weak and the core keeps the spherical collapsed 
shape with the radius i?o- The corona thickness and overall star size increase in this regime upon a decrease in salt 
concentration due to progressively increasing differential osmotic pressure inside the corona. Simultaneously, the 
pulling force applied to the ends of the B-blocks localized at the surface of the core domain increases. At certain 
threshold salt concentration jU(-Ro) becomes equal to the critical minimum force fs, m in given by eq 1321 This point 
corresponds to the onset of the microphase segregation within the core, and appearance of the stretched "legs'. From 
this point on, a decrease in salt concentration leads to a progressive unfolding of the core, an increase in the length of 
the "legs" and, consequently, to an increase in R, D, and H, which is interrupted by abrupt transition from partially 
to completely unfolded core conformation. In the transition point the core size R increases jumpwise and this causes 
a simultaneous jumpwise decrease of the corona thickness H. This drop in the corona thickness occurs because of 
increasing distance between the A-B junction points and decreasing crowding of the A-blocks in the coronal domain. 
However, as a result of these two opposite tendencies the overall (outermost) size of the star D exhibits a small jump 
up. In the law salt regime where the core is completely unfolded, R, D, and H grow again progressively as 0i on 
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FIG. 8: Core radius R (dashed line), corona thickness H (dotted line) and overall star radius D — R+H (solid line) as functions 
of ionic strength. Parameters of the (AB) P star: p — 20, Na = 1000, Nb = 200, va = a 3 /6, \b = 1, C-b = 0.5. Vertical dotted 
lines separate different star regimes. 
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FIG. 9: Diagram of states of amphiphilic (AB) P core-shell star with p = 20, Na = 1000, Nb = 200, va = a 3 /6, \b = 1 on 
{at; 4>ion) plane. ^ and (pl^ are the ionic strengths corresponding to the onset of microphase segregation in the core and to 
complete core unfolding, respectively. 



decreases. 

We see that a progressive decrease in 4>i on causes, in a sequence, two intermolecular conformational transitions: 
The first one corresponds to the continuous onset of the microphase segregation within the core whereas the second 
one is the jumpwise transition corresponding to the complete unfolding of the core (vanishing of the globular phase). 
We denote the salt concentrations at which these transitions occur as and <^'q^, respectively. Their values are 
determined by the parameters of the (AB) p core-shell star (Na, Nb, p, on,, va, Xb), therefore, a multidimensional 
diagram of states can be constructed. Figure|9]and Figure fTOl represent two-dimensional sections of the such a diagram 
of states in coordinates (otb, 4>ion) and (p, (f>i on ), respectively. As follows from the analysis of the first diagram, Figure 
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FIG. 10: Diagram of states of amphiphilic (AB) P core-shell star with Na = 500, Nb = 200, va = 1/6, \b = 1, a% = 0.2 on 
(j>\ <t>ion) plane. and (^Ton are the ionic strengths corresponding to the onset of microphase segregation in the core and to 

complete core unfolding, respectively. 

if the fraction of charged monomers a,b in the corona block A is small (below a certain threshold value ab.min), the 
electrostatic interactions in the corona are not strong enough to induce the unfolding of the collapsed core even at 
vanishing ionic strength of the solution. The value of ah t min is controlled by the hydrophobicity of the core forming 
blocks and can be expressed as 



where fs,mm is given by eg 1321 Moreover, if ctb is only slightly larger than a^minj a decrease in the ionic strength 
stabilizes the microphase-segregated state of the core, whereas the complete core unfolding does not occur even in the 
salt-free solution. 

The diagram of states in (4>i on , p) coordinates, Figure [11)1 allows us to analyze the influence of the number of the 
star arms p on the character of the conformational transitions related to the core unfolding. The diagram presented 
in Figure \W\ shows that <p^ n corresponding to the onset of the microphase segregation in the core is only a weakly 
increasing function of p. By contrast, the characteristic salt concentration corresponding to the jumpwise 

complete core unfolding is a decreasing function of p which approaches the zero at some finite value of p {p « 70 
for our set of parameters). In the stars with larger number of arms the core unfolding occurs continuously upon a 
decrease in the salt concentration. This behavior is related to a decrease of number of monomers in the "minimal 
core" and a decrease in the magnitude of the jump in the A-B junction points position upon transition from the 
microphase-segregated to completely unfolded state as p increases (see Figure \5\ and Section "Core domain"). Hence, 
the microphase-segregated region on the phase diagram in Figure [TU] extends at large p to low salt concentration but 
one should realize that at small 4>i on the core is (almost) completely unfolded. 




(37) 



VI. SCF MODELING 



To check the validity of our theoretical predictions and the correctness of the underlying model assumptions, we use 
Scheutjens-Fleer self-consistent field (SF-SCF) numerical approach to study the structure of the AB core shell star 
and conformational rearrangements upon changes in environmental conditions, i.e., upon change in the ionic strength. 
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A. SCF formalism 



The self-consistent field (SCF) method is well-known for the modeling of inhomogeneous polymer layers. At the 
basis of the approach is a mean-field free energy which is expressed as a functional of the volume fraction profiles 
(normalized concentrations) and self-consistent field potentials. The optimization of this free energy leads for Gaussian 
polymer chains to the Edwards diffusion (ED) differential equation 26 , which for an arbitrary coordinate system may 
be expressed as 

(ks is the Boltzmann constant, T is the absolute temperature, a is the monomer unit size). The Green's function 
G(r, r'; n) used in eq[351 is the statistical weight of a probe chain with the length n having its ends fixed in the 
points r' and r. The self-consistent potential u(r) represents the surrounding of the chain and serves as an external 
field used in the Boltzmann equation to find the statistical weight for each chain conformation. Consequently, Greens 
functions G(r, r'; n) that obey eq|38]are related to the volume fraction profile of the polymer <p(r) or, in the case 
of a multicomponent system, to the set of volume fraction profiles y>fc(r) where k — 1,2, ... specifies the component. 
As long as the potential it(r) is local (i.e., there are no long-range forces) it depends on the local volume fraction 
iys(r). The electrostatic interactions in the system are taken into account on the level of the Poisson-Boltzmann 
approximation: local electrostatic potential is calculated from the Poisson equation, which accounts for the charge 
density distribution created by all charged species (i.e., charged monomer units and mobile ions). This makes up 
the system of self-consistent field equations which is solved iteratively: one assumes an initial volume fraction profile 
(/?(r), then computes the potential u(r), the set of the Greens functions using eg 1381 and derives a new volume fraction 
profile f'{v). The procedure is then repeated until the sequence of approximations converges to a stable solution. To 
solve these equations rigorously, it is necessary to introduce a numerical algorithm. Such numerical scheme invariably 
involves space discretization (i.e., the use of a lattice). Here we follow the method of Scheutjens and Fleer (SF-SCF)2£ 
who used the segment size a as the cell size. A mean-field approximation is applied to a set of lattice sites. This set 
(often called a lattice layer) is referred to with a single coordinate r. The way the sites are organized in layers depends 
on the symmetry in the system and must be preassumed. The approach allows for, e.g., volume fraction gradients 
between these layers. In order to consider an isolated polymer star, it is enough to use a one-gradient version of 
SCF algorithm in the spherical coordinate system, for which r — r. In this case, all volume fraction profiles as well 
as other thermodynamic values depend only on the radial coordinate r, while the mean-field approximation is used 
along angular coordinates. More specifically, within our model p AB diblock copolymer chains are homogeneously 

pinned at the surface of a small sphere with the radius ro = a y / p/(47r) , where \x] denotes is the smallest integer 

not less than x (the so called ceiling function) . AB copolymer nature of the chains is taken into account via the set of 
Flory-Huggins parameters for monomer-solvent interactions, \A, Xb, and monomer- monomer cross-interactions xab- 
The solvent is assumed to be good (athermal) for A monomers, \A — 0, and poor for B monomers, xb = 1- We chose 
Xab — Xb meaning that the incompatibility between A and B monomers is the same as between B monomers and 
solvent molecules. The block A is a positively charged polyelectrolyte, with quenched charge, the fraction of charged 
monomer units is equal to a& (i.e. a& is the fractional charge per A monomer unit). There are two types of salt ions in 
the system: co-ions Na + and counterions CP. SCF calculations were performed using the SFBox program developed 
at Wageningen Uviversity. 



B. Results of SCF modeling 



To make a quantitative comparison with the results of the theory, we have made SCF calculations for (AB) p 
copolymer star with the parameters similar to those used in the previous section (Figures 151 - fl0]> . that is, we consider 
a star with p = 20 AB-diblock copolymer arms (for p — 20 the radius of the grafting sphere is ro = 2), Na = 1000, 
Nb = 200, fraction of charged A-monomers = 0.5, and monomer-solvent interaction parameters xa = (which 
gives for the second virial coefficient va = & 3 /6), and xb = 1. The Flory-Huggins interaction parameters for 
the salt ions are taken the same as for the solvent molecules. The salt volume fraction was varied in the range 
10 _1 < 4>i on < 10~ 4 . The SF-SCF method has obvious advantages, as compared to the analytical theory, in being 
assumption-free: neither strong chain stretching nor equal radial positions of the junction points and terminal monomer 
units are pre-assumed. As a result, the SF-SCF approach provides a more accurate information about density 
distribution in the core and in the corona of the star. On the other hand the use of one-gradient version of the SCF 
approach "smears out" lateral inhomogeneity of the star with microphase-segregated or unfolded core, Figure [3] 
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Since one of our motivations of using the SCF modeling is to check the underlying model assumptions, let us list 
the latters. The analytical theory was developed under assumption that 

1. all hydrophobic B-blocks are equally (but possibly not uniformly) stretched so all the A-B junction points are 
located at equal distance R from the center of the star; 

2. all hydrophilic A-blocks are equally (but not uniformly) stretched so all the A-endpoints are located on the 
equal distance D from the center of the star; 

3. core and corona (correspondingly B and A blocks) are segregated in space: there is a sharp phantom interface 
separating core and coronal domains; 

4. at strong repulsive interaction in the corona microphase segregation in the core occurs; 

5. at strong repulsive interaction in the corona stretched "legs" of the microphase-segregates "octopus" structure 
are uniformly extended; 

6. local electroneutrality approximation for the corona is fulfilled; 

Figure [TT] shows evolution of radial density profiles of A- and B-monomers with decreasing ionic strength. The density 
profile of the swollen corona (A-units density profile, Figure [TTJ a) has a sharp maximum at the interface between 
A and B domains and then decays with increasing r. A decrease in the salt concentration ^ion leads to progressive 
corona swelling: the maximum value of ipA decreases, the profile broadens (i.e. the overall AB star size grows, in 
accordance with Figure [5]) . The density profile of hydrophobic B-monomer units is qualitatively very different from 
its A-counterpart. Figure [TT1 b shows that at high salt concentration (4>i on = 0.1; 0.01) the distribution of the B-units 
inside the core is nearly uniform, with a narrow depletion zone in the layer adjacent to the grafting point and a 
sharp decay at the edge of the core. With a decrease in (f>i on , the pulling force from the coronal A-blocks grows and 
microphase segregation in the core is clearly seen. The profile has two parts: a dense core and a sparse periphery with 
the density much less than that of the core. As <pi on decreases, the globular core becomes smaller and the low-density 
peripheral part extends at the expense of the depleted core. This picture is in accordance with our model assumption 
4 about microphase segregation in the core. 

Another important characteristic is the radial distribution of the number of monomers of the type i, rii(r). It is easily 
obtained from the corresponding radial density distribution ipi (r) by multiplying it by the number of lattice sites in the 
corresponding layer L(r) (which is the volume of the spherical layer) : rii(r) = (pi(r)-L(r) = tpi(r) -47r(ar 2 — a 2 r+a 3 /3). 
Number distribution is the best way for expressing and analyzing the distribution of the end monomers or A-B 
junction points. The evolution of the distribution of A-B junction points, nj(r), is shown in Figure Q2] a. At 
high salt concentration, when the pulling force exerted by the corona is not strong enough to provoke microphase 
segregation in the core, the distribution is rather narrow with a pronounced peak. Our calculations demonstrated 
that the distributions for <fcion= 0.1 and 0.05 (the latter is not shown in Figure [TH a) are indistinguishable, whereas 
the following decrease in (f>i on leads to broadening of the distribution and to progressive shift of the position of the 
maximum towards larger r. By analyzing the shape of the profiles, we can conclude that even though distributions 
are rather broad each profile has one pronounced maximum, and we do not encounter the situation with two maxima 
corresponding to partitioning of the star arms into two groups, i.e., those whose B blocks are completely embedded 
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FIG. 12: Radial distribution of the number of A-B junction points (a) and A-endpoints (b) in AB-core-shell star with p — 20, 
N A = 1000, N B = 200, XA = 0, Xb = 1, a = 0.5 and different values of <t> im . 
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FIG. 13: Radial distribution of the number of B monomer units in AB-core-shell star with p = 20, Na = 1000, Nb = 200, 
Xa — 0, xb = 1, a = 0.5 and different values of 4>i 0n - 



in the core and those whose B-blocks are (partially) expelled from the core and stretched. This finding justifies the 
used by the analytical theory approximation 1 that all the junction points were placed equidistantly from the center 
of the star. 

The distribution of ends of corona blocks, Figure[12]b has one pronounced maximum, with a decrease in <f>i on leading 
to increase of the radius of the A/B interface, R, the distribution becomes asymmetric, biased towards the edge of 
the star. Therefore, the use of assumption 2 is reasonable. Although, as it was noted above, in the one-gradient 
version of the SF-SCF method we cannot see the stretched "legs" , they can be revealed by plotting the radial profiles 
of the number of B monomers per spherical layer. The ng(r) distribution allow us to analyze the local tension in the 
extended part of the core. Uniform "leg" stretching (according to assumption 5) would correspond to a plateau on 
the tib(z) profile. From Figure [13] we see that in our system this situation is better manifested in the case of very 
strong corona stretching force (very low 4>i on )- Deviations from the flat plateau shape of the ns(z) profiles observed in 
Figure [13] (these deviations are also well seen in the plots of Figure [T4]) are due to the finite width of the distribution 
of the A-B junction points. 

Radial distribution of a number of monomer units found at given distance from the center is also an appropriate 
characteristic for estimating the degree of overlap of the core and the corona (hence, we can test the validity of 
assumption 3). These profiles presented together on the same plot are shown in Figured!] (the scale of the plot is 
chosen for better "zoom" of the A/B interfacial region) We see that n>t(r) and ris(r) profiles (shown by black and 
gray lines, respectively) demonstrate a small overlap, the "maximum" value of the ua and tib in the overlap zone 
corresponds to intersection of the profile and is very close to the value of p (in our case p = 20). Hence, the fraction 
of monomers of "guest" A (B) component in the "host" B-core (A-coronal) domain is small, so we can conclude that 
the core and the corona are well segregated, even under the conditions of strong stretching of the inner B-blocks, that 
is the case at low cj)i on values. By summing up the distributions presented in Figure [Ml another interesting feature of 
the star with partially unfolded core is revealed. At strong corona pulling forces (low ionic strength), the distribution 
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FIG. 14: Radial distribution of the number of A- (black curves) and B-monomers (gray curves) in the vicinity of the A/B 
interface in AB-core-shell star with p = 20, Na = 1000, Nb = 200, \a = 0, \b = 1, a = 0.5 and different values of 4>i 0n (shown 
in each panel). 
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FIG. 15: Radial distribution of total number of monomer units in AB-core-shell star with p — 20, Na = 1000, Nb 
Xa = 0, \b = 1, a = 0.5 and different values of 0; on in the vicinity of the A/B interface. 



200, 



of the total number of monomer units, n — ua + Ub, shown in Figure [T5l has a plateau in the vicinity of the A/B 
boundary. This means that the strongly (and homogeneously) stretched "legs" include not solely the parts of B blocks 
extracted from the B-core but also parts of A blocks adjoining the ends of the B blocks. The uniform stretching of the 
arms is typical for starlike polyelectrolyte in the "osmotic" (dominated by the counterions) regime. With an increase 
in r this plateau is followed by the growing part on n(r) profile corresponding to the decrease in the local stretching of 
star arms which is typical for a convex spherical (or cylindrical) polyelectrolyte brush in the salt-dominance regime^. 
At moderate and high ionic strength such a plateau is not observed. 

Finally, to test the validity of the corona electroneutrality assumption 6, local charge density distribution q(r) = 
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FIG. 16: Radial distribution of local charge in AB-core-shell star with p = 20, Na = 1000, Nb = 200, \a = 0, \b — 1, a = 0.5 
and different values of <f>ion- 

<pNa+ ( r ) + a-b^PA^) — if ci- ( r ) was calculated. Figure [TBI shows that in the corona, the local charge density is perfectly 
zero, the distortions of local electroneutrality are observed at the A/B interface only. 

VII. CONCLUSIONS AND OUTLOOK 

In the present paper, conformational transitions in an amphiphilic (AB) p block copolymer star with hydrophobic 
central blocks and polyelectrolyte terminal blocks were studied by means of analytical mean-field theory combined 
with the SF-SCF numerical modeling. The star is in aqueous media and forms under high ionic strength conditions 
a unimolecular micelle with collapsed B-core and extended A corona. Because of repulsive electrostatic interactions 
in the polyelectrolyte corona the core is subjected to the radial extensional force. The strength of the interactions 
in the corona is governed by the degree of ionization of A blocks and by the ionic strength of the solution. From 
the experimental viewpoint, the salt concentration is the most common control parameter. We have shown that a 
decrease in the ionic strength of the solution may provoke the intra-molecular conformational transition related to 
the unfolding of collapsed inner hydrophobic blocks. More specifically, the unimolecular (AB) p micelle undergoes two 
consecutive conformational transitions: The first one is continuous and is associated to the onset of the microphase 
segregation in the B-core. The second one is the complete unfolding of the core which is accompanied by the jumpwise 
increases of the core size and the overall star size and the corresponding jumpwise decrease in the corona thickness. 
Remarkably, the second transition does not occur in the stars with sufficiently large number of arms. We have 
constructed the diagrams of states of (AB) p star in <5>i on — u\, and $i « — V coordinates. The diagrams contain the 
region corresponding to unimolecular micelles with collapsed core, the region corresponding to the partial unfolded 
core-forming blocks and, finally, the region corresponding to the conformation with completely unfolded core and 
extended B blocks. The developed analytical theory is based on two essential approximations: The first one is the 
approximation of local electroneutrality of the corona. Both analytical and numerical Poisson-Boltzmann analysis 
convincingly proved that for reasonable degree of ionization («& > 0.1) the approximation is very accurate for the 
stars with the number of arms p > 10. The second one is the equal arm stretching approximation. The latter is 
well-justified for the unimolecular micelles with collapsed B-core as well as for the conformations with completely 
unfolded core — . The situation, however, is more delicate in the regime with phase segregated (partially unfolded) 
core. In particular, one could expect that there is an analogy to the collapse transition in hydrophobic polyelectrolyte 
stars and brushes: Theories based on the equal stretching approximation^— predicted a jump- wise "unfolding" 
transition. On the other hand, more accurate approache s 33 ! 34 pointed to the co-existence of two populations of chains 
(collapsed and stretched) and smooth variation of the average dimensions in the transition range due to progressive 
re-partitioning of the chains between the collapsed and the extended states. Therefore, we have complemented our 
analytical theory of the unfolding transition in the (AB) p block copolymer stars by numerical calculations made with 
the aid of assumption-free SF-SCF method. The results of the SCF modeling have confirmed the predicted scenario 
of the core unfolding in the amphiphilic (AB) p star and shown that with an increase in electrostatic repulsions in 
the corona, the core unfolding occurs via microphase segregated states. Importantly, the SCF calculations have 
proven that at any ionic strength (i.e., both in the microphase-segregated and and in the unfolded conformations) 
the end segments of the A-blocks and the A-B junction points exhibit unimodal radial distributions, each with a 
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single maximum. This observation proves that the analytical theory based on the equal stretching approximation 
gives a qualitatively correct (artifact-free) description of the conformations and intra-molecular unfolding transition 
in amphiphilic (AB) p star. 

In the present paper we considered (AB) p star copolymers with a quenched fraction of permanently charged 
monomer units in the A-block. Correspondingly, the ionic strength of the solution was chosen as a control parameter, 
variation in which triggers the intra-molecular conformational transitions. In the case of the copolymers comprising 
weakly dissociating (pH-sensitive) polyelectrolyte blocks, the degree of ionization and, as a result, the conformation 
of the copolymer can be affected by variation in pH as well. Moreover, in the case of branched polyelectrolytes the 
degree of ionization of branches is not controlled solely by the value of pH in the buffer, but depends also on the ionic 
strength of the solution and on the polymer conformation. This coupling causes highly non-trivial responsive behavior 
of the pH-sensitive branched (co)polymers. In particular, continuous variation in the ionic strength at pH«pK may 
cause re-entrant folding-unfolding transition. These effects will be addressed in the forthcoming publication. 

Above we have considered very dilute solutions of (AB) p star copolymers which form spherical unimolecular mi- 
celles. At higher concentrations, the aggregation of (AB) P unimolecular micelles with small number of arms may 
lead to formation of intermolecular assemblies. Depending on the lengths of the hydrophobic and polyelectrolyte 
blocks this aggregation may give rise not only to spherical, but also to cylindrical micelles, lamellar or bicontinuous 
structure^. For amphiphilic ionic/hydrophobic diblock copolymers, the morphology of the inter- molecular aggregates 
is controlled by the length of the A and B blocks, the morphological transitions can be triggered by variation in the 
ionic strength or in pH The connectivity of the amphiphilic block copolymers into stars introduces additional topolog- 
ical constraints which may influence the morphologies of the inter-molecular assemblies of (AB) p stars as compared 
to those formed by homologous diblock copolymers ("single arms"). We remark, however, that because of strong 
steric repulsions between polyelectrolyte coronae of the unimolecular micelles, their aggregation may be significantly 
hindered kinetically. Analysis of the inter-molecular assembly of the (AB) p star copolymers is, however, beyond the 
scope of the present work. 
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